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Abstract 

The work ^3] introduces a naive description of dynamics of point vortices 
on a plane in terms of variables of distances and areas which generate Lie- 
Poisson structure. Using this approach a qualitative description of dynamics 
of point vortices on a plane and a sphere is obtained in the works [141 115j . In 
this paper we consider more formal constructions of the general problem of n 
vortices on a plane and a sphere. The developed methods of algebraization 
are also applied to the classical problem of the reduction in the three-body 
problem. 
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1 Classification of the algebra of n vortices on a 
plane 

1.1 Vortex algebra and Lie pencils 

Let us introduce coordinates of vortices r k (x k , Vk), k = 1, . . . , n by complex variables 
z k — x k + iyk, k = 1, . . . , n. The set z = . . . , z n ) defines a vector in a complex 
linear space C n with an Hermitian form 

n 

(z, w) r = ^ TiZiWi, Ti £ R , (1) 
i=i 

where Tj are vorticities. 

The imaginary part of (0) gives the symplectic form corresponding to the Poisson 
bracket {xi, yj} = jr. The Hamiltonian Ti and integrals of motion P, Q, / of vortex 
system can be written as 

1 n 

^ = T^E r ^ ln l^-^l 2 ' ( 2 ) 

n n 

P + «Q = ^2 / T i z i = (z,z ) r , / = ^Til^l 2 = (z,z) r , (3) 

i=l i=l 

here z = (1, 1, . . . , 1) is a constant vector. 

As relative variables we introduce the squares of mutual distances and doubled 
areas of triangles spanned on three vortices (see Fig. 1) as in [TB*] 

Mi m = (T t - rj) 2 = (xi - Xjf + (yi - yj) 2 , 

(4) 

A ijk = (Tj - A (r fe - Ti) = (xj - Xi)(y k - y-) - (x k - Xi)(yj - yi). 
New variables, as it was shown jTHI, generate a linear bracket, which also linearly 
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depends on inverse intensities 1/1^ 



{Mij, M kl } = 4(— S ik - ^-Sjk^Aiji + 4^— 6 U - j^',/)a, ; 
{Mij, A klm } = (—8 ik - (M u - M im + M mj - M jt )+ 



-S u - —S fl ) (M mi - M ik + M kj - M jm )+ 

+ (|T^m - |T^m) (M ki ~ M a + M tj - M jk ), 



{A{j k , Ai mn } — jr(Aj kn — Aj km ) + -j^(Aj k i — Aj kn ) + (5) 

+ T^(Ajfcm — Aj-jy) + -^-(A ifcm — A ifcrl ) + 
1 i 1 j 

_l_ _l_^ ifcn _ Ajjti) + -^-{A ik i — A ikm ) + 

"I - p - (Ajj n Ajj m ) -|- — (Ajj; Ajj n ) -)- — (Ajj m Aiji). 

Let us set invariant relations corresponding to real motions [13 

-FV/fci = Aij k + A ik i — Auj — Aij k = , (6) 

F ijk = (2A l]k f + M% + M% + - 2{M ij M jk + M {j M ik + M jk M ik ) = , (7) 

Relations © mean that the quadrangle spanned on the vortices ijkl can be con- 
structed of triangles by two ways (Fig. 1). Equations (J7J) are Heron formulas ex- 
pressing the area of a triangle via its sides. 

Let us note that for the structure (j5J) the Jacobi identity is fulfilled only on the 
submanifold defined by the first set of geometrical relations (|B|). Below, we consider 
the relations © to be fulfilled on default. Thus the brackets © define Lie-Poisson 
structure, which wee shall call vortex bracket. 

The bracket (j5J) admits the linear Casimir function ^3] 

To define the real type of the corresponding Lie algebra we indicate an explicit 
isomorphism with some Lie pencil jlOj . 

Remark 1. The compactness conditions of the real form |3J), depending on 
intensities, imply that all mutual distances are bounded from above at any moment 
of time. 



3 



We consider a space of skew Hermitian n x n matrices and subspace L in it 
generated by matrices of the form: 



A Zmfc = m 



-1 



(9) 



Mi 



lm 



m 



\ 



—i 



—% 



■ J 



Let us describe the main properties of this subspace which can be verified by a 
straightforward calculation. 
Proposition 1. 

1) The L is a Lie subalgebra, 

2) This subalgebra is isomorphic to u(n — 1), 

3) The center of this subalgebra is generated by the element ^M m ;, 

m,l 

4) The elements Ai ms generate a subalgebra in L isomorphic to so(n — 1), 

5) The following relation holds for any k, I, m, s: 

Afc/ m + Akms + Aiks + A m i s = 0, 

6) Commutation relations in L coincide with the vortex bracket |3J] for the case 
when all intensities are the same (and equal to unity). 

Therefore, Q determines an n-dimensional (unitary) linear representation of Lie 
algebra u(n — 1). There are no such irreducible representations, that is why it splits 
into a sum of the standard n — 1-dimentional representation and one-dimensional 
trivial representation. This splitting is organized as follows. 

Proposition 2. The representation space V = C n splits into a direct sum of 
invariant subspaces V — V\ ® V 2 , where V\ = C n_1 is given by z x + z 2 + • ■ ■ + z n = 
and a one-dimensional subspace Vi = C is spanned on the vector z = (1, 1, . . . , 1). 
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Let us turn to the case of arbitrary intensities Fi, ... , T n . Consider the Lie pencil 
on the algebra of skew Hermitian n x n matrices generated by commutators 



[X,Y] r -i =XT- L Y -YT 
where T _1 is a real diagonal matrix of the form 



-i 



X, 



(10) 



\ 



r 2 



r„/ 



The remarkable fact is that the subalgebra L is closed with respect to the com- 
mutator [-,-]r-i. Thus the family of commutators (fTnj) generates some Lie pencil 
on L. Moreover, restricting the commutator [-, on L we obtain the Lie algebra 
isomorphic to the vortex algebra corresponding to intensities r 1; . . . , r n . Thus 
one can deduce a symmetric isomorphism between the family of vortex algebras and 
a rather simple Lie pencil. 

Using this construction we shall describe properties of vortex algebras. 

Proposition 3. For positive T{ the vortex algebra is isomorphic to u(n — 1). 

Proof. 

The corresponding isomorphism is constructed as follows. From the beginning let 
all intensities be equal to unity. All matrices Q are skew Hermitian and satisfy the 
following property Az = 0, where z is the vector with coordinates (1, 1, ... , 1). 
In other words, this vector is invariant under the vortex algebra. Therefore its 
orthogonal completment, i. e., the hyperplane V± — {z | Y2 z i — 0} * s a ^ so invariant. 

We consider all skew Hermitian matrices possessing this property (that is, map- 
ping z into zero) They form the subalgebra u(n — 1). Thus the vortex algebra in 
case of equal intensities is embedded into u(n — 1). But their dimensions coincide, 
therefore the algebras also coincide. 

This can be shown more explicitly by choosing another basis in space C n . Con- 
sider a basis ei, e2, . . . , e n where ei, . . . , e n _i is an orthonormal basis in the hyper- 
plane V\ (see Proposition 2) and e n = (1/y/n, . . . , 1/y/n). Rewriting matrices 
from the vortex algebra in the new basis one can see that they take the form 



/ 



.4 



n-l 







o\ 



°7 



(12) 



where A n _i is a skew Hermitian (n—1) x (n— 1) matrix. In this basis the isomorphism 
of the vortex algebra with u{n — 1) is evident. A defect of this basis is that this 
basis can't be made symmetric with respect to all vortices. 
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Let us generalize these arguments to the case of arbitrary intensities. Transition 
to the new basis in C n defines the conjugation of matrices in algebra (JSJ) of the form 
A = CA'C~ l , where C is the transformation matrix. (In our case we can assume it 
to be real and orthogonal C^ 1 = C T ). With this substitution the commutator of 
the vortex algebra transforms to the following form: 

AT- l B - BT- l A = CA'C-^CB'C- 1 - CB'C^T^CA'C- 1 = 
= C[A',B'] r ,C-\ 



where 



/ 



1 n-1 



Tin \ 



\7ln • • • 7n-l,n 7n,n J 

Here r / ri _ 1 is the symmetric (n— 1) x (n—1) matrix corresponding to the restriction 
of the form T _1 onto the hyperplane V\. This follows from the relation C^ 1 = C T . 

Taking into account that the last row and column of matrices A and B vanish 
one finds that the vortex commutator transforms to the form 

Let us note that the matrix T' n _ l is positively defined and real, therefore one can 
take the root (r^_ 1 ) 1 / 2 . Then the substitution 



A[ 



n-1 



(K-i)- 1/2 K-i(K-i)- 1/2 



reduces the commutator B' n _^\Y' _ to the standard one. The matrices remain 

skew Hermitian. This argument shows also that the vortex pencil under considera- 
tion is a subpencil of the standard pencil on the space of skew Hermitian matrices. ■ 
The main result following from this construction is that the properties of the 
vortex Lie algebra corresponding to parameters Ti,...,r n are completely deter- 
mined by properties of the bilinear form T' n _ v This form is the restriction of the 
form T" 1 onto the subspace V\ = {zi + • ■ • + = 0} (simply by the signature of 
this restriction). 

Remark 2. Using the method of proof of Proposition 3 it is not difficult to 
show that in general case the vortex algebra is isomorphic to the algebra u(p, q) 
with p + q = n — 1 . 

Let us find conditions for which the algebra L-p is compact, i. e., isomorphic to 
the Lie algebra u(n — 1). The necessary and sufficient condition is that T' n _ 1 is a 
form of fixed sign. In this case there are the following possibilities (we suppose that 
intensities are finite and differ from zero). 
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1) The T 1 is a form of fixed sign, i. e., all are simultaneously positive or 
negative. 

2) The form T _1 has the signature (n — 1, 1), all Tj except for one are positive. 
Indeed the condition of positive definiteness of the restriction on V\ requires 
that the form T" 1 is negative defined onto a one-dimensional orthogonal com- 
plement to the V\ with respect to T This condition is easy to find if one notices 
that the orthogonal complement to V\ in the sense of is a one dimensional 
subspace spanned on the vector t> r = (T 1 , . . . , Y n ). It has the form 

(v r , u r ) r -i = y^^i < 0- 

3) Similarly, in the case of the signature (l,n — 1). 
Finally we have the following 

Proposition 4. The vortex Lie algebra Lr is compact only in the following 
cases: 

1) All intensities have the same sign. 

n 

2) All intensities except for one are positive and ^2 Tj < 0. 

8=1 
n 

3) All intensities except for one are negative and ^2 Tj > 0. 

i=l 

Remark 3. This proposition coincides with the one which was proved in the 
case of three vortices JTffi . 

If the vortex algebra is not "semi-simple", the conditions are defined similarly. 
It happens iff the form T' n _ l is degenerate on V±. From linear algebra it follows that 
this condition equals to the requirement that the orthogonal complement to V\ lies 
in V\. It means that t>r G V\, i. e. 

r x + . . . + r n = o. 

Let us note one more point. The subspace A C L spanned on vectors A [mn is a 
subalgebra for any algebra of the pencil. Its compactness conditions are the same 
as those for the whole L r , i. e., the form r^ l _ 1 should be positive definite. Thus the 
vortex algebra Lr is compact iff the subalgebra A is compact. 
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1.2 Reduction with symmetries and singular orbits 

Now we explain the origin the nature of Lie pencils connected with the vortex algebra 
and describe the (singular) symplectic orbits of the algebra (jSJ) corresponding to real 
motions. Consider a transition to relative variables (@J) from the point of view of 
reduction with symmetries 

Hamiltonian function (0) and equations of motion of vortices are invariant under 
the action of the group E(2). This action can be represented in the form 

g:z^e^z + az , g e E{2) , (13) 

where z = (1, . . . , 1), a = a x + ia y . Parameters <p, a x , a y define the translation 
and rotation corresponding to the element g G E(2). 

The action (JT3jl is the non- Hamiltonian group action jH]. Indeed integrals of 
motion Q corresponding to translations P, Q and to rotations I form a Poisson 
structure that differs from Lie-Poisson structure of the algebra e(2) by a constant, 
cocycle. Obviously this cocycle is unremovable and standard reduction with the 
momentum is impossible (HI 13 122] ■ To execute the reduction in algebraic form we 
use the momentum map in some other way. 

Consider the action of linear transformation preserving the form (•, -)r ©• The 
corresponding linear operators (matrices) form the group which is isomorphic to U (p, q), 
p + q = n and satisfies relations 

xrx + = r, r = dia g (r 1 , r n ). 

Operators of the corresponding Lie algebra A e u(p, q) are defined by 

AT + TA + = . (14) 

Since r + = T from ()14j) . the matrix AT is skew Hermitian. After the substitution 
<f)(A) = AT the standard commutator transforms into the commutator [•, (fTU|) : 

<t>([A,B]) = [</>(A),<l>(B)]r-i. 

The advantage of this approach is that for any Tj the algebra (fT4*j) is represented 
on the same space of skew Hermitian matrices. However, instead of the standard 
commutator it is necessary to use [■, -]r-i- In this case a family of Poisson brack- 
ets {-, -}r-i appears on the corresponding coalgebra u*(p, q). 

Linear vector fields corresponding to operators ((H)) in complex form can be 
written as 

V A = Az = 0(A)r _1 z . 
These fields are Hamiltonian [7] with 

H A = '-(Az, z) r = V^rz = % -(<f>(A)z, z) = V0(A)z . (15) 
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The bracket {•, -jr- 1 of quadratic Hamiltonian functions is 

{^ {A) (z),^ (B) (z)} r - 1 = I([^),0(B)] r - 1 z,z), (16) 

and it corresponds to the commutator [•, -]r-i- 

We can describe the momentum map /i(z) explicitly using the standard identifi- 
cation u(n) and u*(n) with the help of the form (A, B) = TtAB 

Tr/i(z)0(A) = Tr^z T 0(A)z = Tr Qzz T j <p(A). 



Thus 

/ 



% rp % 



u(z) = -zz 
PV ; 2 2 



Z\Z\ ■ ■ ■ Z\Z T 



(17) 



\Z n Z\ • • • z n z n J 

The formula (|17j) defines a map onto a singular symplectic orbit of the algebra u(p, q) 
with the commutator [•, -]r-i- Now the integral (jSJ) / = Yl F Casimir func- 

tion, therefore the reduction with it is carried out. In case of all positive (negative) 
intensities the orbit is topologically homeomorphic to CP™ -1 , because under the 
map ([17]) all points of the form e ltp z (orbits of the rotation group action ([11])) stick 
together. 

Let us carry out the reduction with the remained integrals P, Q ® on the reduced 
space of matrices ([17]) . Due to non-commutativity of P,Q we can reduce only one 
degree of freedom j^J. The constant vector field corresponding to the translation in 
the direction a = a x + ia y is generated by the linear Hamiltonian 

H a = (z,Taz ) = (z,az ) r - 

It's easy to show that Hamiltonians Ha (fTB^) commuting with H a are generated by 
matrices A for which 

ATz = (f)(A)z = . 

So 4>{A) belong to the subspace L defined in the previous section (matrices (fTTj) 
don't belong to this space). 

Now it is easy to see that squares of mutual distances and areas (jH) also admit 
a natural representation of the form 

Mij(x,y) = \zi - Zj\ 2 = i(MijZ,z) , 
A%jk(x,y) = i{A ijk z,z) , 

where My, are matrices ©• The momentum map /x(z) into the corresponding 
algebra u(p', q'), p' + q' — n — 1 has the form 

i ( (z,z )rz \ /_ (z,z )rz 



^=-i-Wn-W • (18) 
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Matrices (fTSjl satisfy the relation /x(z)z = 0, i. e., they belong to the subspace L 
defined in the previous section. In compact case the appropriate orbit is homeo- 
morphic to CP n ~ 2 and corresponds to the reduced phase space (after reducing two 
degrees of freedom). The rank of matrices (j!8|) is equal to the unit. This is their 
characteristic property. 

Remark 4. The structure and compactness conditions of the reduced phase space 
can be found by investigation (by methods of analytical geometry) of a common level 
surface of the first integrals J3J). It is interesting that the algebraic approach also 
allows to solve this pure geometrical problem. 



1.3 Symplectic coordinates 



Let us describe an algorithm of construction of symplectic coordinates for the re- 
duced system of n vortices in the case of compact algebra (u(n— 1)). The method 
of reduction of the vortex algebra (for this case) to the standard representation has 
been described above. In this case elements of the algebra are represented by skew 
Hermitian (n — 1) x (n — 1) matrices with an ordinary matrix commutator: 



/ 



A 



%h\ 



X12 + iyi2 
ih 2 



xi jn _i + iyi,n-i\ 

X2,n-1 + iVl.n-l 



\-£l,n-l + iyi, n -i 



(19) 



/ 



Let us define a Poisson map M 2( - n ^ 
formulas 

h = Pi, 



x 2 ,n-i + W2,n-\ ■■■ ih n -i- 

1) onto a singular orbit u*(n — 1) by 



u in 



h J 



1, n 



1 . 



(20) 



Vij = \fhihj cos(y?i - (pj) , 

The canonical bracket {fi,Pj} = o~ij turns in the standard Lie-Poisson bracket 
on u*(n — 1). The orbit via variables Pi, tpi is written as 

Pi + P2 H 1- Pn-i = C = const . 

This orbit is formed by matrices of rank 1 (2 x 2 minors equal zero). To introduce 
symplectic coordinates we carry out a symplectic transformation 



^0 =Pl +P2 H hPn-1 

Tl =P2 H \-Pn-l, 



50 = l fl, 

51 = ~(pi + >P2 , 



(21) 



r n -2 = Pn-2 , 



Sn-2 — —Pri-l + '■Pn-2 ■ 
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It is easy to see that the map (|2U|) doesn't depend on s , therefore s\, . . . , s„„ 2 , 
r 1; . . . , r„_ 2 define symplectic coordinates of 2(n — 2)-dimensional orbit CP n ~ 2 . 
Remark 5. (Co)algebra (u*)u(n — 1) admits a linear change of variables hi i— > 

n-l 

hi+XD, D = ^2 hi, X = const preserving commutator relations. In this connection, 
i=i 

the orbits given by matrices F e w(ri) of rank 1 map into orbits for which the 
matrix (A — XD) has rank 1. 

The orbit given by Heron relations (in compact case) after the transition from 
variables M, A to standard variables of algebra u(n — 1) turns out among the above 
indicated orbits for some X. 



1.4 Canonical coordinates of a reduced four-vortex system. 
Poincare section 

Using the above described algorithm we shall point out canonical coordinates ex- 
plicitly for a four-vortex problem of equal intensities. 

Consider the coordinates on the algebra u(3) corresponding to the matrix repre- 
sentation of the form (fT^j) 



A, 




%3 + m 

ih 2 
-xi + iyt 




(22) 



In accordance with the above arguments these coordinates are expressed linearly via 
areas and squares of distances 



z), Vi 



E 

k<l 



kl 



■( yi) M kl {z), 



hi 



k<l 



III 



kl 



(23) 



To define coefficients ^.^m^y^m^^ let us use the Proposition 3 and a matrix 
realization (jUJ) of elements A r , M kh where A 1 = A234, A 2 = — A 134 , A 3 = A 12 4, 
A 4 = — A 12 3- We choose the orthogonal matrix C, which reduces matrices of the 
vortex algebra (}2*3*j) by the transformation C 7 AC to the form (fT2*|) as: 



C 



1 

2 
1 

2 



V 



(24) 



Let us set one of coordinates hi,Xi,yi in the obtained matrix A3 is equal to 1 and 
the others are equal to and solve the system of linear equations. Thus we find 
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(25) 



(26) 



corresponding coefficients S^La, m ?4)' m (hi)' ^ n ^ s case 

xi = -(-Ai + A 2 + A 3 - A 4 ), x 2 = i(-Ai + A 2 - A 3 + A 4 ) , 
X3 = -(-Ai - A 2 + A 3 + A 4 ), yi = -{M u - M 23 ) , 
V2 = \ (M 13 - M 24 ) , z/3 = ^ (Mia - M 34 ) , 

/i 4 = i(Mi 2 + M 34 + M 13 + M 24 - M i4 - M 23 ) , 
h 2 = -(M 12 + M 34 - M 13 - M 24 + M 14 + M 23 ) , 
h 3 = -(-M 12 - M 34 + M 13 + M 24 + M14 + M 23 ) . 
Denoting canonical coordinates ()21)1 as (g, G, /i, if), we have 

x x = -y/(H - G)Gsing, x 2 = y/(D - H)Gsin(h + g) , 

X3 = —a/ (£) — H)(H — G) sink, yi = ^ (H - G)G cos g , 
y 2 = ^(D-H)Gcos(h + g), y 3 = y/ (D — H)(H — G) cos h , 
h! = D-H, h 2 = H-G, h = G, 

where D = | ^ w is a constant of the Casimir function «(3), in this connection 
< G < H < D. 

The Hamiltonian in new variables can be represented in the form 

j-{ln({D - G) 2 - A(D -H){H- G) cos 2 h) + 
+ \n((D-H+ G) 2 - 4{D - H)G cos 2 (/i + g)) + (27) 
+ ln(tf 2 -4(H-G) cos 2 ^)} . 

Let us construct the Poincare map on a plane (g, G) for different values of the 
energy E = 7i(h, H, g, G). Surfaces of energy function in figures 2-4 for fixed g = go, 
G — G , f(h,H) = TC(h, H, g,G) show the complex structure of the isoenergetic 
surface. 

We define an intersecting plane by the equation H = H*. The value should 
be chosen so that Vgo, Go < H* equation f(h, H*) = E has the unique solution with 
positive (negative) derivation As is obvious from the presented Fig. 2-4, that it 
is fulfilled not for all (analogous conditions for the intersecting plane h = h* are 
practically never fulfilled). 

The maximal possible energy E? = —4 In 2 ~ — 2.77 . . . corresponds to the Tomp- 
son configuration for D = 1, in this connection H — G — |. In this case the phase 
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portrait on a plane (g, G) consists of the only line G = ~. Phase portraits for the 
energies E < Et are represented in Fig. 5-10. Regions where the motion is impos- 
sible (f(h, H*) = E has no solutions) corresponds to painted regions in Fig. 5-10. 
Figures show that with decreasing of the energy the stochastic layer increases first 
occupying all plane (Fig. 7, 8). Then it decreases and remains only near unstable 
solutions and separatrices (Fig. 9, 10). In the limit E — > — oo one of pairs of vortices 
merges and we obtain an integrable problem that is the system of three vortices. 



1.5 Lax— Heisenberg representation 

As a result of the reduction we can write equations of motion on the orbit of the 
coadjoint representation of the Lie algebra u(p, q). This orbit is singular and consists 
of matrices of the form 



L 



-7*7* 



where Yl z i = ( m a system connected with the center of vorticity). 

According to a general principle (because of the simplicity of the corresponding 
Lie algebra [11J) we can rewrite equations in Lax-Heisenberg form: 



dL 
~dt 



[L,A] 



where A = dH(L) is the differential of the Hamiltonian of the form 

H(L) = Y,^i^sMij(L). 

Here M%j is interpreted as an element of the Lie algebra u(n — 1) and Mij(L) is 
a standard pairing between the algebra and coalgebra. In our case the identification 



2 is fulfilled. 



The explicit form for the differential of the Hamiltonian is given by the formula: 

1 



dH(L) = J2 V ^ 



M tJ (L) M ^ 



i. e., dH(L) is a linear combination of matrices M^. Thus the expression for the 
matrix A has the form 



/V. 



A = dH(L) 



i 1 ! 



r 2 ri 



|z2-Zl| 2 



. r„Ti 



__£lT2_ 

\zi-z 2 \2 

r 2 r j 

\Z2-Zj\ 2 



\zi-z„\ 2 

r 2 r» 

\Z2-Z n \ 2 



"|z„-z n _i| 2 Z-^j \z 2 -Zj\ 2 / 
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If all intensities coincide with each other and equal to 1, the matrix A is simpli- 



fied: 



i 



i 



A = i 



\Z 1 -Z j \' 2 \z 1 -Z 2 \ 2 

i v i 

^3 \*a-zj\ 3 



\Z2-Z X \> 



\z 1 -z n \ 2 

1 

\z 2 -z n \ 2 



\ 



\ \z n -Zi\ 2 



i i / 



In case of different intensities the additional procedure is required to reduce the 
commutator to the standard form. 



1.6 Stationary configurations 



In terms of L — A pair the stationary configuration is naturally described. The 
stationary condition is equivalent to the fact that the matrices L and A commute: 
[L, A] = 0. In our case it means that 



-Azz 1 



-zz T A = --z(Az) T . 

2 2 v ; 



Denoting b = Az, one has bz T = —zb T . It can be shown that it is possible 
only if b = iXz, A G M. It means that z is an eigenvector of a matrix A (with an 
imaginary eigenvalue). 

As a result we obtain a rather natural set of relations 



f^3 1*1-. 



1 



\Z1-Zj\ 2 



\Z2~Z^ 



. 1 

\ \z n ~Z\\ 2 



|*i-*2| 



1 



i |*2-*j| 2 



|*1— *ra| 2 
1 

\Z2-Z n \ 2 



\ 



|2„-2 n _l| 2 |z 2 -Zj| 2 

This can be rewritten in a simpler form 



= Azfc, where k — 1, 



^2 



, n. 



A 



^2 
W 



(28) 



Equations ()28j) can be obtained directly from the equations of motion under the 
condition that every point rotates around the origin with the same velocity A. 

Stationary configurations have been studied in several works (see |TJ [TBI 123] )■ 
Apparently new results can be obtained with the help of the following arguments. 



1. Stationary configurations can be interpreted as eigenvectors of a matrix A. 
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2. Stationary configurations can be interpreted as singular points of Hamiltonian 
on an orbit (which is on CP" -1 ). As a Hamiltonian one can take a function 
H = rifc^j \ z k ~ z j\ 2 - ^ is a positive function which vanishes on a submanifold 
of the "collapse" . 

Using these arguments and by the investigation of a commutativity condition 
[L, A] = 0, it would be interesting to give a theoretical explication of configurations 
from the "Los- Alamos" catalogue. For these configurations stable states of rotation 
realize on concentric circles ("atomic covers" by Kelvin) (see Fig. 11), where for the 
system of 11 vortices two possibilities 11 = 2 + 9 or 11 = 3 + 8) |2] are pointed 
out. As a rule these configurations possess the same type of symmetry (a rotational 
one or have a plane of symmetry). Recently in a short note in "Nature" [3] non- 
symmetric stationary configurations have been indicated for the vortex system of 
equal intensities. 

Theoretically the simplest problem is to determine a number of collinear configu- 
rations in dependence on the ratio of intensities. In celestial mechanics (for positive 
masses rrii) the answer is given by the Multon theorem. In accordance with this 
theorem there is the only (rotating) collinear configuration for any permutation of 
masses mi, . . . , m„ (m* > 0) (the proof of this theorem is contained in j2H])- For 
the three-vortex system the number of collinear configurations depends on the type 
of the algebra defined by Poisson brackets [HJ ITS] . Apparently there is such a 
connection in general case (for the system of n vortices). 



2 Solvable problems of vortex dynamics 

In this section we consider in parallel the problems of dynamics of point vortices on a 
plane and sphere P^HS]. The position of vortices on a sphere in the "absolute" space 
is defined in cartesian r^Xj, y^Zi) or spatial polar coordinates (9i,<fi, % = 1..N) 

Xi = R sin 9i cos , y\ = R sin 9i sin ipi , Zi = R cos 9i , 

where R is the radius of a sphere. Let us remind (see J3|), that the Hamiltonian 
function and Poisson bracket are defined by relations 

H = -— Vr^lnLRW^, 

%<3 

{ipi,cos9j} = ^|T%> i,j = h---,N, 

where 9^ is the angle between z-th and j-th vortices, and Tj is a vortex intensity. For 
relative variables on a sphere the same labels — M^, are used, however, their 
meaning is different, than for a plane JH]- So My is a length of a chord connecting 
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the vortices, and is equal to the ratio of volume of a tetrahedron spanned on 
three vortices and the center of a sphere to ^R. 

Mij =(xi - Xjf + (yi - yj) 2 + (zi - Zj) 2 , 

A 1 f ^ (29) 

2.1 The particular case of the problem of n vortices, the 
reduction to the problem of (n — 1) vortices 

There is a special case of the problem of n vortices on a plane and a sphere, for that 
case the problem can be reduced to the system of (n — 1) vortices with the same 
algebra of Poisson brackets, but reduced Hamiltonian function. The procedure of a 
reduction in Hamiltonian exposition corresponds to restriction of the system on a 
Poisson submanifold [TH], which is defined in this case by involution conditions of 
integrals of motion. For vortices on a plane the integrals are given by relations (J5jl. 
necessary conditions of reduction accept the form 

N 

^r i = 0, Q = P = 0. (30) 
1=1 

For a sphere integrals of motion are obtained similarly: 

N N N 

F\ — R Ti sin 6 1 , cos cpi , F2 = R Tj sin 0i sin cpi, F3 = R Tj cos Q; L , 
i=i 

and conditions of reduction 



i=i i=i i=i 

(31) 



F 1 = F 2 = F 3 = . (32) 

Remark 6. The equations H3U\) have the following geometrical meaning: each vor- 
tex is situated at the center of vorticity of all remaining vortices. Really, expressing 

N-1 

for example Tn = — T,, from we have found that the absolute coordinates 

i=i 

of the N—th vortex are defined by expressions 

N-1 

1=1 



N-1 

i=i 
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The geometrical meaning of relations \3ty on a sphere is a little different, it means, 
that the center of vorticity of a system of n vortices coincides with the geometrical 
center of a sphere. 

The mentioned below reduction happens because there is a possibility of calcula- 
tion of positions of n vortices basing on positions of (n — 1) vortices. 

For a determination of additional invariant relations for relative variables we 
shall use identities j2H| 

1 - 

- T k(M jk - M ik ) = P{ Xi - Xj ) + Q( Vi - Vj ) (33) 
k=i 

N 

for plane which are fulfilled only under condition of ^^r« = 0, that is also valid for 
the case of sphere jTHj *=i 



- T fc (M ifc - M lk ) = R(F 1 (x i - x,) + F 2 ( Vi - Vj ) + F 3 ( Zi - Zj )) . (34) 
k=i 

Using conditions ([3*0J). we have obtained Cf^ invariant relations for 

N 

J2r k (M jk -M ik ) = 0. (35) 



k=l 



These relations we shall supplement by relations for A^, which also follow from 
(|3TIj). ((nil)- It is possible to show by a straightforward calculation that the complete 
set of these relations determines a Poisson submanifold with the structure (j3J) and 
its generalization for a sphere j!3j . 

Using representation for Mij in absolute coordinates (jlj) and (j^Jj) . it is simple 
to verify that among the equations ()35|) only n — 1 are linearly independent. Using 
these equations, we can express squares of distances from one of the vortices (nth 
vortex) to all the others M kn , k — 1, . . . , n — 1 via mutual distances between n — 1 
vortices M^, i,j = 1, ...,n— 1. Then substituting these expressions into initial 
Hamiltonian, we shall obtain a system with a bracket of the problem of n — 1 
vortices and reduced Hamiltonian function. 

For the case of four vortices the explicit formulas of squares of distances from 
the first three vortices up to the fourth vortex M 14 , M 24 , M 34 have a form 

. . r|M 13 + t\m 12 + r 2 r 3 (M 12 - m 23 + m 13 ) 
" Mu = (r. + r 2 + r 3 r ■ (36) 

the reduced Hamiltonian is obtained by substitution (}3l)|) into initial Hamiltonian 

(ED- 



IT 



The system of three vortices is integrable (independently on a Hamiltonian given 
on algebra (j5J)), the described case corresponds to a special case of integrability of 
problem of four vortices (originating in work by Kirchhoff [21] )• The phase portraits 
of this system at various intensities are presented in the previous work [T5] . 

Remark 7. The indicated special cases of integrability correspond to a situation 
at which one of integrals reaches the extreme value. It is obvious, that in this case 
the system necessarily has additional invariant relations. For integrable systems that 
situation produces an additional degeneration. Delone case for Kovalevskaya top is 
one of examples. In this case the Kovalevskaya integral, which is the sum of two full 
squares will vanish, and the two-dimensional tori degenerate into one- dimensional 
torus (periodic and asymptotic solutions). 

2.2 Particular solutions in the problem of four vortices 

General equations of motion of vortices with some restrictions on intensities Tj 
allow a finite group of symmetries, the elements of the group are the permutations 
and reflections in some planes. Such discrete symmetries do not produce the general 
integrals of motion and do not allow reducing the order of a system. However, 
presence of these symmetries is connected to invariant submanifolds. The solution 
on that submanifolds can, as a rule, be obtained in quadratures 

Let us consider two problems of dynamics of four vortices on a plane and a sphere 
possessing two types of symmetry — central symmetry and reflection symmetry (for 
a plane, central symmetry and axial symmetry). 

a. A centrally symmetric solution at D = 0. Equations of motion of four 
vortices on sphere (flat case is obtained with the passage to the limit R — > oo) with 
condition r\ = r 3 , T 2 = T 4 allow invariant relations 

M u = M 23 = M u = Mi , M 12 = M 34 = M 3 , ^ 
A234 = A124 = Ai , A213 = A314 = A2 . 

(The relations (J37|) don't define the Poisson submanifold) . 

The equations ()37|) have the following geometrical meaning: a centrally symmet- 
ric configuration of vortices (the parallelogram), keeps this symmetry in all moments 
(see Fig. 12). 

The analysis of a centrally symmetric solution in absolute variables through 
explicit quadratures is carried out by D. N. Goryachev [T7] (see also |3J). However, 
he has not made clear the qualitative properties of motion and has introduced the 
very confusing classification. Let us introduce the qualitative analysis of relative 
motion. 

The equations describing evolution of sides Mi,M 3 and diagonals M 13 = M 2 , 
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M 24 = M 4 of the parallelogram, have a form 



Mi 


= 4r x A 2 1 


^M2 


1 

M 3 


M 


= 4r x A 2 1 




1 

M 2 


M 2 


= 8r 2 A 2 1 




1 

Mi 


M 4 


= srxAi 1 


V^3 


1 

Mi 



U L. 

^M 4 M 3 



(3? 



For a sphere (for a plane i? — > 00) the geometrical relations between M, A can 
be written as 

Ai (AR 2 - M 2 ) + A 2 (AR 2 - M 1 - M 3 ) = , 
Ai (4i? 2 - Mi - M 3 ) + A 2 (4i? 2 - M 4 ) = . 

System is solvable with respect to A 1; A 2 under condition of 



(39) 



2(M l + M 3 ) - (M 2 + M 4 ) + ^- (M 2 M 4 - (M 1 + M 3 ) 2 ) = . (40) 

The linear integral of the equations (jHBj) corresponding to the Casimir function (jHJ) 
has a form 

d = 2r 1 r 2 (M 1 + m 3 ) + r 2 M 2 + r 2 M 4 . (4i) 

With the help of relations ()39p and regularization of time, the system (138(1 can 
be reduced to two inhomogeneous equations describing evolution of sides of the 
parallelogram Mi, M 3 . 

For simplicity we consider the limited case D — 0, that is also a necessary 
condition of a collapse [24J. Geometrical interpretations on a plane and a sphere are 
a little different, therefore we shall consider these cases separately. 

1) Plane. 

Relations (|3T? |) .(f4T) |) for a plane (R — > 00) have a form 

A 1 = -A 2 = A, 

42 

2(M X + M 3 ) - (M 2 + M 4 ) = 0. 

Taking into account, that in (j4"Tj) D = 0, we found the equation, describing a 
trajectory of system on plane Mi, M 3 (sides of a parallelogram) 

rfMi Mi (2r 1 r 2 (M 1 + m 3 ) + (r 2 + r 2 )M 3 ) 

dM 3 m 3 (2r 1 r 2 (M 1 + m 3 ) + (r 2 + r 2 )Mx) ' 1 j 
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The solution of this equation has a form 



(M 1 - M 3 ) 2 = (Mx + M 3 ) 2 -C(M 1 + M 3 ) 
a = ^ + p 2 , C = const . 

41 ll 2 7 



(44) 



(Index a < 0, because with £> = intensities r 1; T 2 have different signs). 

In a quadrant M 1 > 0, M 2 > (which corresponds to physical area) three types 
of trajectories (JUj) are possible depending on a: 

1°. —1 < a < - - all trajectories are closed, started from the origin of co- 
ordinates, tangenting axes OMi OM 2 (see. Fig. 13, a), a) (inhomogeneous 



2°. a < — 1 — all trajectory are curves asymptotically approaching to coordinate 
axes (see Fig. 13, b). 

3°. a = —1 - - all trajectories are straight lines moving through the origin of 
coordinates (see Fig. 13, c) (homogeneous collapse) 

The physical area on a plane Mi, M 3 is defined by the part of a positive quad- 
rant (Mi > 0, M 3 > 0), for which the inequality A 2 > is fulfilled. When the 
trajectory reached a boundary A 2 = 0, the sign in the equation should be 
changed (reflection) that corresponds to the same trajectory passing in the opposite 
direction. It is easy to show, that the equation A 2 = is defined by two straight 
lines on a plane Mi, M 3 , which are situated inside a quadrant Mi > 0, M 3 > for 
any intensities T\, r 2 . Hence, except for a case 3° vortices are moving in the bounded 
area without collisions and scattering. In the system in case 3° the homogeneous 
collapse of all vortices or homogeneous scattering happen. 

For the problem on a plane each point of the trajectory on the plane Mi,M 3 
correspond to two configurations of vortices, which differ only by permutation: 1 <-> 
3, or 2 <-> 4 (see Fig. 12). 

Remark 8. The condition of a homogeneous collapse (a = — 1) from the analy- 
sis of motion in absolute variables is obtained in \2$ . As result of quasi-homogeneity 
of the equations the condition of a homogeneous collapse can be obtained by investi- 
gation of existence conditions of solutions of the form Mi = Crf^, G\ = const of the 
system. 

Remark 9. For a plane the system h38(l with a regularization 



collapse). 



dr 



4A 



d,t 



M 1 M 2 M 3 M 4 
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can be reduced to a homogeneous Hamiltonian system 

= r 1 M 1 M 4 (M 3 - M 2 ) + T 2 M 1 M 2 (M 4 - M 3 ) , (45) 

dr 

= T 1 M Z M A {M 2 - M x ) + T 2 M 2 M 3 (M 1 - M 3 ) , (46) 

dr 

= 2T 2 M 2 M 4 (M 1 - M 3 ) , (47) 

dr 

^ = 2T 1 M 2 M 4 (M 3 - Mi) (48) 
dr 

with a Poisson bracket of the form 

{M u M 2 } = l-MM.M.M, , {Ma, M 3 } = ^ ^ - i-J M X M 2 M 3 M± , 
{M 2 , M 3 } = -^M 1 M 2 M 3 M 4 , {Mi, M 4 } = -^-MiM 2 M 3 M 4 , (49) 

1 1 1 2 

{M 2 , M 4 } = 0, {M 3 , M 4 } = —M 1 M 2 M 3 M± 



r 



2 



and Hamiltonian 



h = 2r!r 2 in Mi + r? in m 2 + 2rir 2 In M 3 + In M 4 . 

The rank of a bracket is equal to two, hence dividing the bracket on MiM 2 M 3 M 4; 
it can be reduced to a constant without violation of the Jacobi identity. 

2) Sphere. 

For a centrally symmetric configuration on a sphere we also consider a projection 
of trajectories to a plane Mi, M 3 (see Fig. 14). The form of physical area defined by 
inequalities Af > 0, % = 1, 2 varies. Moreover, there are two differences from a case 
of plane connected with the nonlinearity of the equation (|40p. At first, each point 
Mi, M 3 correspond to two solutions of the equation (pUj) . and therefore two various 
(spatial) parallelograms with the given sides on a sphere, these parallelograms are 
not reduced to each other by permutation 1 3 (2 «-> 4) . Secondly, the system (JHS|) 
is not reduced to two quasi-homogeneous equations. 

Let us express Ai and A 2 from (JH9*j) with preservation of homogeneity 

Ai = (M 4 - Mi - M 3 )A , A 2 = (M 2 - M x - M 3 ) A . (50) 
Also we shall make change of the time 

dt = - - , , dr . 

MiM 3 

So we obtained the system describing evolution of variables Mi, M 2 , M 3 , M 4 
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_ . . f r {M 2 -M 1 -M 3 ){M 3 -M 2 ) {M A -M l -M 3 ){M 3 -M A ) \ 

]l - U| V ~M~ 2 • ' 2 M A ) ' 

„ L (M 2 -M 1 -M 3 )(M 2 -M 1 ) (M,-M 1 -M 3 )(M i -M 1 ) \ 

M 3- M 3^ — + r 2 — j , 

M' 2 = 2r 2 (M 2 - Mi - M 3 )(M! - M 3 ) , 
M4 = 2r!(M 4 - Mi - M 3 )(M! - M 3 ) , 

(51) 

here primes denote derivation on new time r. 

Let us introduce a linear change of variables reducing system (J51|) to two differ- 
ential equations, which has the most simple form. 

Let us choose variables x, y, M, N: 

x = T.M2 + T 2 M 4 , y = T\M 2 - r 2 M 4 , 

M = M x + M 3 , N = Mi - M 3 . ! ? 

With the help of a relation (|41|) where D = and (J4*U|) which have a form 

y 2 

x = "i6i? 2 r 1 r 2 (r 1 + r 2 )' (53) 

we exclude x, M from (|51|). In outcome we obtain the equations of evolution y, N 
4(T 1 + T 2 )yN, 

1 r 2%-r 1 (r 1 + r 2 ) z,^^ 
" Wi + r 2 ) 3 r?r| U^-r^r. + r,)) l 41 ^ 1 + r *) * " 

- y 2 (2% - T 2 + T 2 2 )(2ky(T 1 + 2T 2 ) - T 1 (T 1 + 3r 2 )(ri + T 2 ))) + 

+ 2% + r2( / ri + r2 \, (41*1*^ + r 2 ) 4 iv 2 - (54) 
r 1 (% + r 2 (r 1 + r 2 ))V 1 21 2j 

- y\2ky - Y\ + T 2 )(2%(r 2 + 2I\) + T 2 (r 2 + 3^)^ + r 2 ))) 

hereinafter k = 1/R 2 . 

The projection of trajectory y(r), N(t) on the plane M 1; M 3 is defined by the 
formulas M x = |(M + N), M 3 = |(M - JV), where 

r r 2 
M ^ 2r 1 r 2 2 (r 1 + r 2) y+ i6gr 1 r 2 (r 1 + r 2 p - ( r ' r * <0 >- < 55 > 

The equation (jo3j) allows real solutions only for M < M ma:r = — ^r-j^ R 2 - There- 
fore, in case of a sphere the physically accessible area on a plane is defined by 



dy 

dr 
dN 

dr 
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inequalities 



M 1 + M 3 < M max , 

2(M 1 M 4 + M 3 M 4 + M X M. A ) - Ml - Ml - Ml - ^-M X M 3 M± > , 

R 

2{M 1 M 2 + M 3 M 2 + M X M 3 ) - M 2 - M 3 2 - Mf - -^-MiM 3 M 2 > . 



(56) 



4A? 



(57) 



4A 2 



In the area (J57|) the equation (J55|) has two roots, therefore we have on a plane 
M\, M 3 a projection of two different areas of possible positions of vortices defined 
by inequalities (Joljj) and equations (jB3j) . If the areas (JoTlj) do not reach a straight 
line Mi + M 3 = M maa; (Fig. 14 b, c), the point, beginning to move inside one area, 
remains there in all moments. When the point reached boundaries (A, = 0), it 
is necessary to change a sign of time, and the trajectory is passed in the opposite 
direction. In a case when the straight line Mi + M 3 = M max passes through the 
interior of areas (JSTIj) (Fig. 14 a) the trajectory, that had reached that line, passes 
from one area in another and is described by other solution of the equation (jB3j) . 
Characteristics of phase portraits in case of a sphere is the occurrence new motionless 
points, absent in a flat case. If the second necessary condition of a collapse is realized 
(r 2 + T 2 , = — 4r!Fi), we can see on a Fig. 14 c), that some of trajectories starting 
at the origin of coordinates goes to origin again and do not reach the boundary of 
area (collinear positions) and some trajectories goes back to origin only after the 
reaching of boundary. It means, that the collapse remains homogeneous only near 
the origin of coordinates. 

b. A reflective-symmetric solution. For a system of two cooperating vortex 
pairs Ti = — 1^, T 2 = — T 3 the general system of four vortices on a sphere and on a 
plane also allows invariant relations 



These equations have the following geometrical meaning. The vortices located at 
the initial momentum in tops of a trapezoid will form a trapezoid during all time of 
motion (see Fig. 15) The given configuration has a reflection (axial) symmetry. 
The equations describing the evolution of sides and diagonals of a trapezoid have a 




Mi 3 = M 24 = M 3 , 

A124 = Ai 3 4 = A 2 . 



(58) 
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form 
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1 


W A 




1 


1 


W A 


~ m7 


1 
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Mi 




1 


1 


Mi 


"M3 



M 3 = 4r 2 Ai (4r ~ -vr ) + 4r 2 A 2 
M 2 = 8r 2 A 2 



1 


1 


M~ 2 


~ Ms 


1 


1 


Ma 


~ Mi 



(59) 



M 4 = 8riAi 

Here M 2 = M i4 , M 4 = M 23 — of the bases of a trapezoid. 

The geometrical relations between M, A in this case have the identical form on 
a plane and sphere: 

AiM 4 + A 2 (M 3 -M 1 ) = 0, 

Ai(M3-M!) + A 2 M 2 = 0, [ ' 

the condition of their resolvability with respect to Ai, A 2 has a form 

M 2 M 4 - (M x - M 3 ) 2 = . (61) 

The integral of the momentum has a form 

d = 2rir 2 (Mi - m 3 ) - r^M 4 - r^M 2 . (62) 



From the equations (J59H6TJ) it follows that trajectories of a system in space M 1? M 2 , M 3 , M 4 
coincide for cases of a plane and sphere. The difference between these problems con- 
sists in form of physical areas defined by inequalities (J56j) . 

With the help of the equations (j6Ti|) we shall express A 1; A 2 by 

A, = (M 2 - Mi + M 3 )A, A 2 = -(M 4 - M x + M 3 )A , 

and also we shall exclude A with the help of regularizing substitution of time 

4A , 
dt = , r dr . 
M 1 M 3 

To reduce the regularized equation to the system of two equations, we make the 
substitution (|52j). As a result of the existence of an integral of the moment between 
the variables x, y the following relation is fulfilled 

4r x r 2 „ 2(rir 2 ) „ . . 

Dx + y 2 +-y-^Dy + D 2 = 0. (63) 



Ti + r 2 Ti + r 2 

Using (jpHj) and equation (JBTJ, we find 



y 2 D 

N = y - + . (64) 

^t 1 t 2 d 4rir 2 v ; 
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With D ^ from relations and (ffii jl all M 1? M 2 , M 3 , M 4 can be expressed 

through M, y. For the variables M, y we shall get the equations of the form 



dy fa 2 -D 2 )fa(ri+r 2 )+D(ri-r 2 )) 

dr 4riT 2 D ' 

+8(D-y) 2 DrlM-(D 2 -y 2 ) 2 ]- (65) 

- 32^1^ + y) ^ " Fl) " y(r ' + r M 16 ™ M ' + 

+8(D + yfDY\M - (D 2 - y 2 ) 2 ] . 

By analogy with a centrally symmetric solution, we shall consider a projection of a 
trajectory to a plane M 1; M 3 . Besides inequalities (J56*|) the physical area is defined 
by an additional condition (corollary of the relation ([64))) 

£>(£> -41^2^) > 0. (66) 

Two various points of a phase space (two solutions of the equation (|51jl) are projected 
in each point on a plane M 1; M 3 , satisfying to inequalities (|56|). It is more 

evident if we imagine that two different areas defined by an inequality (|56j) are glued 
together along a line Mi — M 3 = N . When a point reached a boundary Aj = 0, it 
is reflected back and moves on the same trajectory, and at reaching boundary (|66|1 
the point passes from one area to another (see. Fig. 16 a, 16 b). In figures for 
convenience we have unrolled the areas, which have been glued together along a 
line Mi — M 3 = N , where N is a solution of the equation D — 4r 1 r 2 A r = 0. The 
difference of a plane from a sphere is exhibited in form of physical areas. For a plane 
(see Fig. 16 a) two types of trajectories exist: 

1. Trajectories are touched once boundaries Aj = and then leaving to infinity; 

2. Trajectories stayed between boundaries Aj = 0. 

Different motions of vortices correspond to these cases: in the first case the vortices 
run up, passing only once through a collinear configuration, in second case vortex 
pairs alternately pass one through another (leapfrog of Helmholtz) , remaining on 
bounded distance. Only trajectories of the second type exists for a sphere. The 
analysis of leapfrog on a plane was carried out by other method in [23]. 

Let us introduce also, for completeness, the equation of motion and geometrical 
interpretation at the zero moment D = (necessary condition of a collapse [21] )■ 
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According to (|63|) in this case y = 0, and all mutual distances Mi, . . . , M 4 can 
be expressed via variables x, M by formulas 

M — i\/f _ Fix i\T _ x 

m z - r 1 (r 1 +r 2 ) > m * ~ r a (ri+r a )> iy ~ r 1 +r 2 > 
Afi = i(M + JV), M 3 = |(M — AT) . 

The equations of motion for m = (I\ + r 2 )M and x have a form 
dx 



-2x 2 



dr 

dm i (v x r 2 \ 2 i m r 2 \ 2 

— — = m — 2xm x . 

dr 2 Vr 2 rj 2 Vr 2 rj 

The trajectory, defined by the system (|87|). is obtained from equation 



(67) 



g^ 4/a - 1 r 2 \ 
m = x ^A73i' a= 2U + r\J' (68) 

where C is a constant of an integration. The equation, defining the form of area (Aj > 
0) on a plane m, x has a form 

'< ro - m >- 8g(r 1 + r,)0 >O - <69> 

Analyzing (|68|). (|69|) near to the origin of coordinates m = x = it is possible to 
conclude that for reflectional symmetrical solution the simultaneous collapse of four 
vortices is impossible. 

The discussed above integrable systems of vortex dynamics allow the complete 
enough analysis with the help of qualitative research of dynamical systems on a 
plane. At the same time application of classical methods of an explicit solution with 
the help of the theory of special (Abelian) functions produces very bulky expressions, 
and that expressions are not permitting to make any of representation about the 
real motion J7| . 

3 Algebraization and reduction of the three-body 
problem 

The reduction in the three-body problem of celestial mechanics is for the first time 
systematically considered in the lectures of Jacobi [HI]. He dwelt upon barycentric 
coordinates, which permit to ignore uniform motion of the centre of masses, and also 
to the procedure of elimination of momentum (elimination of a knot). Construc- 
tively this process reducing in a plane (spatial) case in a system with three (four) 
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degrees of freedom was executed by Radout, Brunce, Charlier, Lie, Levi-Chevita and 
Whittaker, the review of their researches contains in the |27| . The most natural 
and complete reduction of order in the three-body problem belongs to van Kampen 
(E. P. van Kampen) and Wintner (A. Wintner) [20]. This problem is also considered 
in the work [21] • All these classical results are based on the theory of canonical 
transformations and connected with bulky calculations. Here we give more geomet- 
rical method of the reduction of order in a flat three-body problem, which could be 
easily generalized to a spatial n-body problem. It is connected with a preliminary 
algebraization of the reduced system and subsequent introduction of canonical coor- 
dinates on a symplectic orbit. In contrast with the classical approach this procedure 
uses (algebraic) Poisson structure of the reduced system only and has no connec- 
tion with Hamiltonian (and therefore it can be derived for other potentials, which 
depend on mutual distance between bodies). The algorithm of reduction under our 
consideration is universal enough and uses only algebraic methods. Undoubdedly, 
the study of the algebraic form of n-body problem (not mentioning more natural 
symmetrization) has large prospects for celestial mechanics. 

3.1 An algebraization of a system 

The Hamiltonian of a plane three-body problem can be written in the form 

*=5E*H2><i r «- r 'i> (70) 

t=i i^j 

where Ti(xi,yi), Pi(pi X ,Piy) are two-dimensional vectors of position and momentums 
of particles, for which components (pi X , Xi, pi y , yi) of Poisson brackets are canonical. 

Let us choose new variables, which characerize relative dynamics of particles as 
the quadratic forms by canonical variables: 

1. Quadrates of mutual distances 



Mi = 


l r 3 


r I 2 
— r 2| , 




M 2 = 


I 1 *! 


-r 3 | 2 , 


(71) 


M 3 = 


l r 2 


-r x | 2 . 





2. Scalar products of momentums and neighbouring sides of a triangle formed by 
points (see Fig. 17) 



x 2 = 


(ri 


- r 3 ,Pi), 


A 3 = 


(r 2 


- ri,pi) 


Y x = 


(r 3 


- r 2 ,p 2 ), 


Y 3 = 


(r 2 


- ri,p 2 ) 


Zi = 


(r 3 


- r 2 ,p 3 ), 


z 2 = 


(ri 


- r 3 ,P 3 ) 
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Since the Hamiltonian (J7U|) is dependent only on absolute quantity of momen- 
tums and mutual distances, it can be expressed through relative variables (f7TH72|) . 
Moreover these variables form a Lie-Poisson bracket 



{X 2 ,X 3 } 


— x 2 + x 3 , 


{^2,^1} 


= 0, 


{X 2 ,Y 3 } = 


-Y x - Y 3 , 


{x 2 ,z 1 } 


= X 2 + x 3l 


{x 2 , z 2 } 


= —X2 — Z 2 , 


{X 3 ,Y 1 } = 


—X2 — X 3 , 


{X 3 , Y 3 } 


= X 3 + Y 3 , 


{x 3 , z x } 


= 0, 


{X 3 ,Z 2 } = 


Zi + z 2 , 


{Yi,Y s } 


= -Y x - Y 3 , 


{Yx,Z x } 


= Y x + Z 1: 


{Yi,Z 2 } = 


-Y x -Y 3 , 


{Y^Z,} 


— —Z\ — Z 2 , 


{Y 3 ,Z 2 } 


= 0, 


{Zt,Z 2 } = 


Zt + Z 2 , 


{X 2 ,M X } 


= 0, 


{X 2 ,M 2 } 


= -2M 2 , 


{X 2 ,M 3 } = 


M 1 -M 2 -M 3 , 


{X 3 ,M 1 } 


= 0, 


{X 3 ,M 2 } 


= -M x + M 2 + M 3 , 


{X 3 ,M 3 } = 


2M 3 , 


{Y 1 ,M 1 } 


= 2M l5 


{Y!,M 2 } 


= 0, 


{YuM 3 } = 


-M 2 + M 1 + M 3 


{Y 3 ,M 1 } 


= M 2 -M 1 - M 3 , 


{Y 3 ,M 2 } 


= 0, 


{Y 3 ,M 3 } = 


-2M 3 , 




= -2M l5 


{Z U M 2 } 


= M 3 -M 1 - M>. 


{Z U M 3 } = 


0, 


{Z 2 ,M 1 } 


= -M 3 + Mi + M 2 , 


{Z 2 ,M 2 } 


= 2M 2 , 


{Z 2 ,M 3 } = 


0. 


{M U M 2 } 


= 0, 


{M 3 ,M 2 } 


= 0, 


{M U M 3 } = 


0. 

(73) 



It can be shown, that variables M, X, Y, Z commute with an integral of the total 
angular momentum of a system with respect to an arbitrary point and quadrate of 
total momentum. Hence, mapping (|71H72|1 corresponds to the reduction with these 
integrals, while the rank of an initial canonical bracket (equal 12) reduces four units 
down. Thus rank of a bracket (J73j) is equal to eight, and the Casimir function is the 
quadrate of total momentum (which is expressed through relative variables (j71fJT2*j) 
in contrast to the angular momentum). 

Remark 10. The method of the algebraization of a dynamical system which was 
used here and in the vortex dynamics is a special case of a general method. It is 
based on the fact that the quadratic functions of variables p, x form the algebra spin) 
with respect to canonical symplectic structure u = dpi A dxi There exists 

a certain subalgebra in sp(n), and its generators commute with integrals of motion 
of a system to each particular problem (Hamiltonian) . 

3.2 Barycentric coordinates and Poisson submanifolds 

The above mentioned reduction concerns an arbitrary inertial system, for which two 
projections of total momentum and total angular momentum form in general case a 
noncommutative set of integrals. In the system of the centre of inertia (barycentric 
frame) the total momentum is equal to zero, hence set of integrals is commutative 
and the reduction for one more degree of freedom is possible. To reduce the order 
in the algebraic form it is necessary to restrict a system on a submanifold of zero 
total momentum in the algebra (|73j) . 
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We choose new generators of the algebra J73J). They correspond to eigenvectors 
of Killing form jH] 



Si = ^= (X 2 + X 3 + Y 1 + Y 3 + Z 1 + Z 2 ) 



S3 

s 5 



1 



(2X3 - Y 3 + 2Z 1 -Z 2 + X 2 - Y x ) 



4^3 

X 2 — Yi — Y 3 + Z 2 , 



Ni = -(M x + M 2 + M 3 ) , 

y6 



No 



V2 



S 2 

(Mi - M 3 ) 



(Y 3 + Z 2 -X 2 -Y 1 ) 

(X 2 — X 3 — Yi + Y 3 

-X 2 — X 3 + Yi + Z\ , 
1 



Z x - 



N 3 



(Mi - 2M 2 + 
(74) 



The variables S$, Sq are proportional to projections of total momentum on two 
sides of a triangle (see Fig. 17). A linear span S$, Sq forms an ideal of algebra (|73j) . 
hence submanifold of zero momentum 



S, = 



Sa = 



(75) 



is Poisson submanifold. Restricting a system on (|75jl and ordering remaining vari- 
ables as follows x = (S±, S 2 , S3, S4, Ni, N 2 , N 3 ), we obtain the table of commuta- 
tors 








S3 


-S2 








iV 3 


-N 2 \ 




-S3 





-Si 





-iV 3 









S2 


Si 








N 2 


Ni 



















-Ni 


-N 2 


-iV 3 







N 3 


-N 2 


iVi 













-N 3 





-JVx 


iV 2 













N 2 


iVi 





iV 3 








/ 



(76) 



The Casimir function of structure ()76|) coincides with quadrate of the total angular 
momentum with respect to the centre of masses 



Ml 



(S,N) 2 
(N,N) 



(77) 



where (a, b) = a^bi — a 2 b 2 — a 3 b 3 is a scalar product in the Minkowski space. 
Numerator and denominator of this fraction are Casimir functions of a subalge- 
bra so(l, 2) ® s R 3 with generators (Si, S 2 , S3, Ni, N 2 , N 3 ). Moreover (N, N) = 2A 2 , 
where A is the square of a triangle formed by particles. 

The algebra l-j, corresponding to a bracket ()76|) is the semidirect sum of subal- 
gebra formed by elements Si, i = 1, ... ,4 and three-dimensional commutative ideal: 
h = (so(l,2)@R) @ s R 3 . 



29 



The quadrates of momentum and mutual distances can be written as 



Pk 



1 (e fc , N) gj - 2 (e fc , S, N) £ 4 + 2 (e fc , S) (S, N) (e k , N) (S, S) 
3 (N, N) 

M fc = (e fc ,N), fc= 1, 2, 3, 



(78) 



where S = (Si, S%, S3), N = (Ni, N 2 , N 3 ) and vectors have the form 




2,-y/3,l), e 2 



V6 



1 



(-2,0,-2), e 3 



1 



2,-^3,1) 



and satisfy the relation (ei,ej) = 1 — 5ij. Here (a, b, c) is a determinant matrix 
which is formed by components of vectors a, b, c. 
The Hamiltonian has the form 



3.3 Orbits and symplectic coordinates 

In algebra (|7fi|) the orbits appropriate to physical motions of a system are divided 
in two types. 

1 . Regular six-dimentional orbits are surfaces of a level of Casimir function (J77|) . 
The quadrate of area of a triangle is nonnegative, therefore for physical sym- 
plectic sheets (N, N) > 0. Topologically such orbit is diffeomorphic TL 2 x 
]R + x R. Here TL 2 is a tangent bundle of a pseudosphere (N, N) = ci, 
its "radius" has non-negative values c x 6 IR + = [0, 00), and the last factor 
corresponds to a linear span S±. 

2. It is possible to show (analyzing a rank of a matrix ()76|)). that four-dimentional 
orbits diffeomorphic to tangent bundle to a cone (N, N) = pass through 
points satisfying the equations (N, N) = 0, (S, N) = 0. These orbits corre- 
spond to the motion of particles along a straight line (the area A is equal to 
zero). 

Remark 11. The orbits passing through points with negative value of quadrate 
of area (N, N) < or with zero quadrates of distances have no physical sense. 

We construct symplectic coordinates on a regular orbit similar to Anduier-Depri 
coordinates in the rigid body dynamics 12.. For this purpose we will apply the 
following algorithm. 

For algebra ()7f)j1 we shall consider the following consequence of the inserted 
subalgebras so(l, 2) C so(l, 2) © s R 3 C I7. 



H 




(79) 
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On a subalgebra so(l, 2) we shall choose L = Si as an action. For a Hamiltonian 
field of vectors with the Hamilton function 7i = L 



dSi _ dS 2 _ dS 3 _ 

dl ' dl ' dl 2 

time parameter along an integral curve Z is a desired angular variable canonically 
conjugate L: {l,L} = 1. Choosing constants of an integration so that the relation 
{£>2, S 3 } = —Si can be fulfilled, we obtain 



Si = L, S 2 = VL 2 -G 2 cosl, S 3 = VL 2 -G 2 sinl. (80) 



We shall choose the Casimir function G = y~S\ — S\ — S 2 of subalgebra so(l, 2) 
as a Hamiltonian 7i = G on a subalgebra so(l, 2) © s R 3 . Then we integrate the 
equation that is linear with respect to iVj 

dNi _ S 3 N 2 - S 2 N 3 dN 2 _ S 3 Ni - S X N 3 dN 3 -S2N1 + S X N 2 
dg G dg G dg G 

We find dependence of constants of integration on L, I, G from commutators of the 
subalgebra so(l, 2) (B s R 3 - 



N^^ + y^J^-A^cosg 



N 2 = ^ VL 2 - G 2 cos / + k J% - A 2 cos g cos I - J % - A 2 sin 3 sin / , (81) 



jV 3 = ^L 2 - G 2 sin Z + §y fl - A 2 cos # sin Z + y g| - A 2 sin # cos Z , 

where if = (S, N) , A = (N, N). 

We shall choose as a last action variable S — S4. As a result we obtain the 
following expressions 



. M Z L V£ 2 - G 2 M 2 , 



^2 = e s (-^jjV L 2 -G 2 cos Z + ^f- - lcosgcosZ - y ^ - lsingsinZ), ^ 

iV 3 = e s ( VL 2 - G 2 sin Z + | - 1 cos g sin Z + ~ 1 sin # cos z) , 

$4. s . 

here M z is a constant of total momentum with respect to a system of the centre of 
masses. The variables S%,S 2 , S 3 express by the formulas (|50|) . 

It is possible to introduce the other pair of canonical variables x, y, instead 
of s, S, by the formulas e s = x, S = xy. 
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Remark 12. Using the above mentioned algebraic structure of three-body prob- 
lem, it is also easy to reconstruct the results of classics in a problem of reduction 
of order. However, obtained expressions will remain bulky enough. Remind that in 
researches of Brunce, Whittaker and van Kampen as positional variables the mutual 
distances are chosen. Reduction of the order which has been carried out by Charlier, 
uses distances of three bodies from the general centre of inertia. 

Using the algebraic form of the equations of the three-body problem it is easy to 
investigate its particular solutions. One of such solutions was founded by Lagrange. 
In this case all of three bodies are situated in tops of an equilateral triangle (which 
in an even more special case does not change sizes). The second solution belongs 
to Euler and defines collinear stationary configurations. All these solutions can be 
defined by a method, explained in |2E1 12H1- and also one can obtain a system of 
appropriate invariant relations which are not determined, Poisson manifolds. 

It is easy to show (see EH]), that there are only two "rigid-body" stationary 
configurations in the three-body problem — Euler and Lagrangian. Collinear config- 
urations are investigated in detail in the n-body problem. By the Multon theorem 
their number is equal to % |2o] . 

The topological proof of the Multon theorem is given with the use of the Morse 
theory. In the paper [26J series of hypotheses about the amount of noncollinear 
configurations in a n-body problem are also expressed. As far as we know the 
majority of them are not proved until now. We hope that the algebraic form of the 
reduced system offered in this application, will allow to achieve further progress in 
this problem. 

The autors are grateful to N. N. Simakov for the help in numerical computa- 
tions. In Russian outcomes of this paper are introduced in the book A. V. Borisov, 
I. S. Mamaev "Poisson structures and Lie algebras in Hamiltonian mechanics". 
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Figure 2: G = 0.3, g = tt/2 




Figure 4: G = 0.3, g = 
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Figure 12: 





b) 



0.2 0.4 0.6 Mj 
ri = -0.183, 

r 2 = 2.o. 



c) 



0.2 0.4 M 
T 1 = -0.183, 
T 2 = 0.683. 



Figure 13: 
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Figure 16: 
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Figure 17: 
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